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ABSTRACT 

We present time-series echelle spectra of the Pa/? line of the T Tauri star SU Aur, ob- 
served over three consecutive nights. The line shows strong variability (~ 10 per cent) over 
the velocity range (100 km s^^, 420 km s^^) in the red broad absorption component, and 
weaker variability (^ 2 per cent) over the velocity range (—200 km s^^, km s^^) in the 
blue wing. The variability in the velocity range (—200 km s^^, km s^^) is correlated with 
that in (200 km s^^ , 400 km s^^), and the variability in these velocity ranges anti-correlates 
with that in (0 km s~^, 100 km s^^). The mean spectrum from the second night shows a sug- 
gestion of a blue-shifted absorption component at about —150 km s^^, similar to that found 
in the Ha and H/3 lines. We find the position of the subpeak in the red absorption component 
changes steadily with time, and its motion modulates on half the rotational period. We also 
find that the modulation of the line equivalent width is associated with a half and a third of the 
rotational period, which is consistent with the surface Doppler images of SU Aur Radiative 
transfer models of a rotationally modulated Pa/3 line, produced in the shock-heated magneto- 
spheric accretion flow, are also presented. Models with a magnetic dipole offset reproduce the 
overall characteristics of the observed line variability, including the line equivalent width and 
the motion of the subpeak in the red absorption trough. 

Key words: stars:formation - stars: individual: SU Aur - circumstellar matter - stars: pre- 
main-sequence 



1 INTRODUCTION 

' T Tauri stars (TT S) young (<~ 3 x 10® yrs, 

lAppenzeller & Mundll Il989^ low-mass stars, and are thought 
to be the progenitors of solar-type stars. Classical T Tauri stars 
(CTTS) exhibit strong Ha emission, and typically have spectral 
types of F-K. Some of the most active CTTS show emission in 
higher order Balmer lines and metal lines (e.g. Call, H and K). 
They also exhibit an excess amount of continuum flux in the 
ultraviolet (UV) and infrared (IR). 

Line profile studies of CTTS show evidence for both 
outflows and inflows, as seen in the blue- shifted absorption 
features in Ha profiles (e. g. iHerbid Il962h and the inverse 
PCygni (IPC) profiles (e.g. iKenvon et alJll994l: lEdwards et"ai] 

Il994h respectively. Typical mass-loss rates of CTTS are about 
lO"'^ Mp, yr~^ to 10"^ M p, yr~^ (e.g. lKuhilll964 lEdwards et alj 

Il987t iHartigan et alj[T99^ . and mass -accretion rates range from 
10"^ Mp, yr~^ to 10"^ Mp, yr~^ (e.g.lKenvon & HartmanrJ 19871: 

iBertout. Basri. & BouvieJl988HGullbring et ajl998h . 

In the currently favoured model of accretion on to 
CTTS, the accretion discs are disrupted by the stellar mag- 
netosphere which ch annel the gas from the disc on to the 
steUar surface (e.g. lUchida & Shibatj Il985t iKoenigll ll99lL 
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ICoUier Cameron & Campbelll Il993t IShu et alj Il994ii) . This 
picture of the accretion flows is supported by the obser- 
vation that CTTS have relatively strong (~ 10^ G) mag- 
netic fields (e.g. [J ohns-Krull, Valenti. Hatzes. & Kanaad 
1 19991 lOuenther & Emerson .1996: .Svmington et aL .2004^ 
and by radiative transfer models which reproduce the 
gross characteristics of the observed profiles for some TTS 
(e.g. Muzerolle, Calvet, & Hartmann 2001). 

SU Aurigae is a bright (J ~ 7, IChakrabortv & Gell2004 
young star with a G2 spectral type iHerbia [l960|), a nd has 
been classified eith er as a CTTS (e.g. lOiampapa et alj Il993t 
lBouviereta'i]|l993h or early-type TTS (e.g.'Herbst et al."'l994j)- It 
is also the prototype 'SU Aur type' iHerbig & Be ll 1988). Based 
on modelling of the spectral energy distribution (SED), SU Aur is 
known to have an accretion disc (~ 400 au) along with an outer 
envelope (e.g. 
observations of 



Akeson et al]|2002h. Near-infrar ed c oronagraphic 



Nakaiima & Golimow ski (1995^ and lGradvet"ai] 



y001|) have shown that the object is associated with a refle ction 
nebula and outflows. Most recentlv. lchakrabortv & Gji2004b have 
presented high resolution (0.3 arcsec) J-, H- and if-band adap- 
tive optics images displaying nebulosity with a disc-like structure 
around SU Aur. 

SU Aur has been a su bject of a numbe r of va ri ability studies, 
mainly in the optical, e.g. Giampapa et alj il993h . ljohns & Basril 
jl995l) . |Petrov et al.Hl996h . loiiveira et al.H200(]li and most recently 
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IUnmhetalJj2004 . The Ha line, which is most commonly used for 
variability studies of TTS, is produced over a large circumstellar 
volume, and shows evidence of outflow (winds) as well as mag- 
netospheric inflow, [johns & Basri ( 1995) studied the variability of 
Ha and H/3 profiles of SU Aur, and found that the modulation of 
the wind and the accretion components in the profiles are approx- 
imately 180° out of phase. This led them to suggest that the mag- 
netic field axis is slightly tilted with respect to the rotation axis of 
the star. 

The topology of the magnetic field is still in dispute. Despite 
of the overall success of the dipole field geometry used in ex- 
plaini ng many observed features of CTTS, Johns-KruU & G af ford 
y002) re-examined the validity of this assumption. They found 
that observations did not support the relationships between mass, 
radius, and rotation period that magnetospheric accretion theory 
predicts under the assumption of a dipolar field, but found bet- 
ter agreement when a non-dipole field topology was used. In ad- 
dition, for SU Aur. fUnruh et alJ 1200 4) found a lack of the circu- 
lar polarisation across some strong photospheric lines; once again 
in conflict with the expectations for a large-scale dipole. Further, 
lOliveira et alj i200(i) found the time-lagged behaviour in the vari- 
ability of some optical lines, and proposed that this can be caused 
by the presence of an azimuthal component in the magnetosphere. 

Near-IR profiles, such as Pa/3 and Br7 are likely to become 
important probes of accretion in CTTS and their lower-mass coun- 
terparts, since they appear to suffer from less contamination by out- 
flo ws than Ha (e.g. |F olha 1 998: Folha & Emerson 2 001. although 
see IWhelan et alj|2004 and are less affected by dust obscuration, 
enabling the study of accretion in highly embedded objects. Here 
we present time-series Pa/3 spectra of SU Aur, which we use to 
probe the accretion geometry. In Section |2| we give details on the 
observations and the data reduction. The results and analysis are 
presented in Section|3| Radiative transfer models are computed in 
an attempt to explain the overall observed variability feature in Sec- 
tion|4| The conclusions are given in Section|5| 



2 OBSERVATIONS 

A total of 503 Pa/3 (A = 1.28181 /im) spectra with a signal-to- 
noise ratio (S/N) of ~ 90 per pixel (in the continuum) were ob- 
tained over 3 nights (2002 December 1-3) at the United Kingdom 
Infrared Telescope (UKIRT) on Mauna Kea, Hawaii. The approxi- 
mate temporal sampling rate was one spectrum every 3 minutes. 

The Cooled Grating Spectrometer No. 4 (CGS4), a 1-5 jira 
multi-purpose 2-D grating spectrometer with a 256 x 256 InSb 
amy, with the echelle grating and long camera (300 mm), was 
used. This provides a velocity resolution (at 1.28181 /im) of ^ 
16 kms~^ using ~ 1 pix slit width without the Nyquist stepping. 
The wavelength range of the echelle spectra is about 3600 kms^^ 
providing enough continuum on either side of the line for accu- 
rate rectification. Unfortunately there are insufficient photospheric 
features to perform a veiling correction (but see Section l431 . Tel- 
luric standard stars were observed throughout each night (every 2-3 
hours) to correct for instrument response and atmospheric transmis- 
sion. 

The data obtained with CGS4 were reduced using a standard 
procedure (c.f. J'uxlev et a l. 1993): 1. b ias subtraction, 2. flat-field, 
and 3. optimal extraction faornelll98d) . All spectra were reduced 
in the same manner, using the ORAC-DR pipeline program. Intrin- 
sic absorption components in the standard stars, used in step 3 to 
correct for the telluric lines, were subtracted from the original spec- 



tra by using the least square fit of the spectra with a combination 
of two Gaussian functions and a third order polynomial. After the 
preliminary reduction by the pipeline was finished, a wavelength 
calibration was performed using a combination of OH sky emis- 
sion lines and comparison arc lines. Finally we applied heliocentric 
velocity corrections and shifted the spectra into the rest-frame of 
S U Aur, adopting the r adial velocity measurement (-|-16kms~^) 
of lHerbig&Beirjl988l) . 



3 RESULTS 

A general overview of the dataset is presented in Fig. Q The 
mean spectrum shows a classic inverse P Cygni (IPC) profile, 
with the blue-wing of the emission component extending to 
—200 km s^^ and the absorption component with maximum depth 
at +100 km s^^ and extending out to +200 km s^^. The peak in- 
tensity of the line, relative to the continuum, is about 1.25. Al- 
though c omparable in strength, th e Pa/3 profiles of SU Aur pre- 
sented bv lFolha & EmersonllioOll) show significant emission red- 
wards of the absorption dip, and also have maximum absorption oc- 
curring at Okms~^ i.e. they appear to be shifted by ^ 90kms~^ 
blueward compared to our dataset. 

The mean spectra from each night are shown in Fig. |2| It is 
immediately apparent that the blue side of the profile is more stable 
than the red, with variability of less than 10 per cent (in contin- 
uum units) at the line peak and blueward. There is a suggestion of 
blue-shifted absorption at about —150 kms~^ on the second night. 
Interestingly, the b lue-shifted absorption feature at —150 km 
is also seen in Ha jGiampapa etaljl993ft and H/3 jjohns & Basril 
1 199^ profiles, and the modulation period of the absorption feature 
is identified as the rotational period of SU Aur (--^ 3 d). The red 
side of the mean profiles shows strong variability, extending red- 
ward to over +400 km s~ ^ . On the final night there is emission on 
the red side of the absorption trough, and the profile more closely 
resembles that presented by Folha & Emerson ( 2001). 

The greyscale dynamic spectra (Fig. reveal the gradual 
changes that occur over a timescale of hours. On the first night 
there is a feature that accelerates monotonically redward across the 
absorption trough. On the second night both the red and blue sides 
of the profile are lower than the mean. Finally, on the third night 
we see the growth of an emission feature in the absorption dip. The 
characteristic level of this red-wing variability is about 10 per cent. 



3.1 Temporal variance spectra 

In order to quantify the level and significance of the variability seen 
in the Pa /3 profiles, t he temporal variance spectrum (T VS) analysis 
method proposed bv lFullerton. Gies. & Bol'torljl99d) was applied 
to the time series spectra from each night separately, and also to 
all three nights combined. This method statistically compares the 
deviations in line features with those of the adjacent continuum. In 
other words, the level of deviation at a given wavelength is com- 
puted with a weight function which is inversely proportional to the 
signal- to-noise (S/N) level o f the continuum. By using the defini- 
tion of lFullerton et alj \l99^ . the j-th velocity bin of the TVS, in 
the limit of high photon-count, can be written as: 

JV 

1=1 
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Figure 1. The time-series spectra of Pa/3. The mean spectrum (normahsed) 
of the 503 spectra obtained during the three night (2002 December 1-3) is 
shown at the bottom. In the middle panels, the quotient spectra (divided by 
the mean spectrum) as a function of time (HJD — 2452600) are plotted 
as greyscale images with the colour scaled from 1.1 (white) to 0.9 (black). 
The temporal deviation spectrum (TVS)^/^ from the 503 spectra is shown 
on the top panel (c.f. Section ITT] . In all plots, the horizontal axes are in 
velocity space in the stellar rest frame. The greyscale images clearly show 
variability on both red and blue sides of the spectra, but little variability is 
seen near the line centre. The dashed line in the top panel shows the sta- 
tistical significance at the 1 per cent level, and it indicates that the velocity 
range (—200 km s^^, 420 km s~^) has significant variability. 



where Sij is the signal of the j-th velocity bin belonging to the i-th 
spectrum of a time-series, and A*" is the number of observed spectra. 
Sj is the weighted mean spectrum defined as 



1 



Sj — — y ^ WiSij 
1=1 



The weighting factor Wi for the i-th spectrum is defined as 



(2) 



(3) 



where aic and akc are the noise in the continuum for the i-th and 

fc-th spectra. 

As mentioned by iFullerton et alj il996h . a more convenient 
quantity to be plotted is the square root of TVS (temporal devia- 
tion spectrum) because (TVS)^'^^ scales linearly with the spectral 




Velocity (km s ) 

Figure 2. Mean spectra from each night: The first night (solid), the second 
night (dashed), and the third night (dotted). While most of the variability 
occurs in the red side (50 km s^^, 500 km s~^) of the line, small yet 
noticeable changes in flux level can be seen in (-200 km s-l , km s'^). 
The flux level at ~ 250 kms~^ is below the continuum on the first and the 
second night, but it is in emission (above the continuum) on the third night. 



A hint of a blue-shifted absorption component at ' 
the second night. 



150 km s ^ is seen on 



deviations; hence, it provides a more accurate impression of the 
deviation from the variability with respect to the continuum. 

Using equation Q the temporal deviation spectrum of the 
times-series spectra from each night was computed, and the results 
were placed in the top panels of Fig. In the same figure, the 
spectra divided by the mean spectrum of each night (the quotient 
spectra) and the mean spectrum are shown as a function of time. 
Similarly, (TVS)^''^ computed from all three nights is placed in 
the top of Fig. Q 

The (TVS)^''^ value which corresponds to one per cent sta- 
tistical significance (the x^-probability) is also shown in the same 
figure as a reference point. In the absence of real variability, the 
observed (TVS)^'''^ values should be below this level for 99 per 
cent of time. From the (TVS)^'''^ plot in Fig. Q we confirm 
that the line displays strong variability within the velocity range 
(100kms"\ 420 kms"^) in the red wing, and weaker variabil- 
ity within the velocity range (—200 km , km s~^) in the blue 
wing. The (TVS)^''^ plots in Fig.|3|show that the velocity ranges 
in which the variability is most active change from night to night. 

Using the mass ( M, = 2.25 Miti) and the radius (i?» — 
3.6 Rq) of SU Aur from lCohen & Kuh?il979h . we find its escape 
velocity to be V^sc = (2GM*/ii*)^^^ = 490 kms"^ This value 
is an upper-limit for the velocity of the material as it impacts on 
the stellar surface, and indeed we see that the extent of the red ab- 
sorption in the line profile is ~ 420 km ^ , and that the maximum 
velocity of significant variability is ~ 420kms~^. 



3.2 Motion of a subpeak in tlie absorption trough 

In the red absorption trough of the mean spectra (see the bottom 
panels of Fig.|3}, we find the presence of a subpeak at V =~ 250, 
~ 180 and ~ 260 km s^^ for the first, second and third nights. 
We also find the position of the subpeak (indicated by arrows in the 
same figure) changes during one night. It appears that a main cause 
of the variability seen in the red-wing (as seen in the greyscale im- 
ages in the same figures) is due to the motion of this subpeak. 
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Figure 3. The summary of observations from each night (the first night to the third night from left to right respectively). For each night, the nonnalised mean 
spectrum (bottom), the quotient spectra (divided by the mean spectrum) (middle) with colour scaled from 1.08 (white) to 0.92 (black) as a function of time 
(HJD — 2452600), and the temporal deviation spectrum TVS~^/'^ (top) in units of 1/100 of the continuum flux (VoFc) are shown. The dashed line in the top 
plot indicates the statistical significance for 1 per cent level. In the bottom plots, the positions of the subpeak in the red absorption troughs are indicated with 
aiTows, and they seem to be moving from night to night (see Section lXll . 



The positions of the subpeak were estimated for each spec- 
trum, by applying a quadratic fit to the data points around the sub- 
peak. The results are shown in Fig. |4| as a function of time. The 
amplitude of the change in the subpeak position is relatively large 
on the second night compared to those of the first and the third 
nights. To demonstrate the motion of the subpeak occurred dur- 
ing the second night, we plot the time-averaged (At — 70 min) 
spectra in Fig.|5] The figure shows a weak subpeak moving from 
V ~ 170km s"^ to ~ 220km s"^ during the course of the night. 

The subpeak positions were fitted with a function (/) in the 
form of the Fourier series, 

f(^) = Y+ E {^') + {^') } 

n-1 

where t and P are the time of observation and the rotational period 
of SU Aur respectively, and Timax = 3. The reasons for choos- 
ing the form of the fitting function above are as follows: 1. a sim- 
plest possible form is desirable, and 2. the first harmonic terms in 
the Fourier series might be present if the variation is related with 
the stellar rotation, 3. the second harmonic terms in the Fourier 
series might be present if the variation is cau sed by the tilted- 
axis magnetospheric accretion model (e.g. see IShu et alj 119943 : 
Ijohns & B asri^l995l). and 4. the third harmonic terms in the Fourier 
series mig ht be present because the surface Do ppl er images of 
SU Au r bv lPetrov et'ai]<199d) . lMennessied (l993) and lUnruh et"ai] 
j2003) weakly suggest the presence of cool spots on the surface at 
three different longitudes which are approximately equally spaced. 

The fit with rimax = 3 is shown in Fig.|4]along with the fits 
with rimax ~ 1 and rimax ~ 2 for com parison. In all cases, the 
fixed period of 2.7 d jUnruh et alj2004 was adopted. See TableQ 
for the summary of the fitting parameters. The figure clearly shows 
that the fit with rimax = 1 fails to represent the data. On the other 
hand, the fits with rimax = 2 and rimax = 3 are significantly better. 



The two lines (rimax = 2 and rimax = 3) are very similar 
to each other as we can see from the figure and the values of the 
Fourier coefficients. The contributions of the third harmonic terms 
(as and bs) are very small compared to the first and the second 
harmonic terms. To test whether the third harmonic terms should 
be included in the fi tting functions, we have performed an F-test 
(c.f lBevingtonlll969h . From the table, the numbers of degree of 
freedom are V2 = 497 and = 495 for rimax = 2 and rimax = 3 
respectively. Using the reduced chi-square values in the table, we 
find J'^ = {x'('W = 2)-x^(nmax = 3)}/xS, = 2 and 
F ~ X^2/x^3 = 1- Since > F, the third harmonic terms 
should be included in the fitting although their contribution is rela- 
tively small. 

From this analysis we found that the motion of the subpeak 
in the red absorption trough is possibly associated not only with 
P, but also with P/2. This is consistent with the the tilted-axis 
magnetospheric accretion model mentioned earlier. 

In spite of the excellent fits of the data with Fourier series, we 
should not ignore a possibility that the rather large amplitude of the 
change in the subpeak position seen on the second night was caused 
by a single episodic event (e.g. by a temporal and local enhance- 
ment of the accretion flow). Assuming conservation of energy, the 
free-fall time (tg) of an object from the co-rotation radius (r — Rc) 
of the star to the stellar surface (r = R,) can be written as 

where Ksc {Rc) = (2GA/./i?c)^''^, and R, and M, are the stel- 
lar radius and mass respectively. By evaluating the integral, one 
obtains 

.«.^^{|+.^/^(l-.)^/^-arcsin(0} (6) 
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Figure 4. The positions of tlie subpeak (circles) in the red wing of the Pa/3 
profile from the three nights are plotted as a function of time. The data were 
fitted with the Fourier series (equation with the terms up to rimax = 1 
(dash-dot line), nmax = 2 (dashed line) and rimax = 3 (solid line). The 
con'esponding fitting coefficients are summarised in Table [Tl It is clearly 
shown that the line with nmax = 1 does not fit the data very well, but the 
lines with nmax = 2 and nmax = 3 fit the data equally well. The contri- 
bution from the third harmonic terms in the Fouiier series is veiy small. The 
motion of the subpeak in the red absorption trough is associated not only 
with the rotational period of SU Aur (P), but also with half of the period 
(P/2). A typical size of the uncertainty in the peak positions is indicated on 
the lower-left corner. 



where q = R^/Rc- Using the mass ( M^ = 2.25 Mp) and the ra- 
dius (J?. = 3.6 Rq) of SU Aur from lCohen & Kuhii fl979l) along 
with the co-rotation radius (Rc = 3 i?, ) in equation |6| we find 
tff = 10.5 h. This time-scale is comparable to the observing time 
span of a given night, so the association of the subpeak motion on 
the second night with a single episodic event cannot be ruled out 
by timescale arguments alone. However, the systematic trends in 
the acceleration of the feature on the first and third nights strongly 
suggests that the feature is produced by rotational modulation. To 
exclude the possibility, the object must be observed for a few rota- 
tional periods. 



3.3 Auto-correlation map 

The auto-correlation map for the Pa/3 spectra (with all 503 spectra) 
was calculated to examine if variability at a given velocity bin is 
correlated with those at a different part of spectra. The result was 
placed in Fig. |S|as a greyscale image. The correlation coefficient 
value idj) of the i-th row and the j-th column in the map was 
computed by using: 
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where Smn is the signal of the n-th velocity bin in the m-th spec- 
trum of the time-series, and A'^ is the total number of spectra. Sn 
is the mean signal of the ?i-th velocity bin in the time-series. The 
map is symmetric about the diagonal (i = j pixels). The range of 
the correlation values varies from —1 to 1 which correspond to a 
strong anti-correlation and a strong correlation respectively. 
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Figure 5. The time-series spectra of the Pa/3 from the second night (solid) 
and the quadratic fits for the subpeak positions (dashed). Each spectrum is 
the average of 27 consecutive spectra (At = 70 min) in chronological order 
from the top to bottom. Each spectrum is vertically displaced by -0. 1 from 
the previous spectrum for clarity. The position of a weak subpeak in the 
red-wing appears to be moving from V ~ ITOltm 8"^ to ~ 2201<m s^^ 
during the course of the second night. The positions of the subpeak were 
estimated using the oiiginal spectra (before being averaged over 27 spectra), 
by applying a quadratic fit to the data points around the subpeak (Fig.l4l. 
The arrows indicate the measured positions of the subpeak. 



The map is useful to visually identify veloci ty bin regions 
that correlate or anti- correlate with each other (e.g. ljohns & Basra 
1199.4 lOliveiraetalJEoOOi) . The figure shows the profile vari- 
ability in the velocity range (— 200 kms~^, Okms~^) weakly 
correlates with that for (200 kms~^, 400 kms~^). On the other 
hand, they seem to weakly anti-correlate with the variability in 
(0 km , 100 km s~^). The pattern seen in the ma p is very simi- 
lar to t hat seen in the H/3 auto-correlation function of lOliveira et alj 
i200Cl) . They also found three similar wavelength ranges in which 
the flux levels conelate and anti-correlate with one another. 

The auto-correlation maps with time lags of P/2, P/3 and 
P/4 where P i s the rotational pe riod were also computed (with 
P = 2.7 d from lUnmh et aj2004 . Because of the poor time cov- 
erages (3 nights at a single location) of the data, we were unable to 
draw any significant conclusion from those maps, and they are not 
shown here. 



3.4 Line equivalent width 

The line equivalent widths (EWs) of Pa/? profiles (shown in Fig.0 
are computed by averaging/rebinning over 10 profiles obtained 
consecutively (approximately over 30 minutes), in order to de- 
crease the size of variance due to the relatively low average S/N 
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Table 1. The summary of the fitting parameters in equation|3]for the positions of the subpeak shown in Fig.|2] A'^ and xt ^re 
the number of data points and the reduced chi-square respectively. 
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Figure 6. Auto-correlation map of the Pa/3 spectra (3 nights) with no time 
delay. The stronger the correlation, the brighter the pixels in the map. Sim- 
ilarly, the the stronger the anti-con'elation, the darker the pixels are. The 
flux changes seen between (-200 km km s^^) and (200 km 
400 km s~^) are correlated with each other. On the other hand, the flux 
changes seen between those wavelength ranges weakly anti-correlate with 
the flux changes seen in (0 kms-i, 100 kms-l). 



(~ 90) in the original profiles. Three different ranges of velocity 
bins were used for computing the EWs. The total EWs were com- 
puted using the velocity range -500 km s"^ <V < 500kms"\ 
The EWs of the red wing (the red EWs) were computed using 
km s^^ < V < 500 km s~^, and those of the blue wing (the 
blue EWs) with -500 kms"^ < V < kms"^ The results 
are shown in Fig. During the first night, the total EW changes 
between —0.1 and —0.8 A. It changes between +0.1 and —0.4 A 
on the second night, and during the third night it changes between 
—0.7 and —1.9 A. The negative sign in the EW values indicates 
that the line is in emission. From night to night, the average EW 
changes from —0.4 to —0.2 A, and then to —1.5 A. 

As we can see from this figure and also from the temporal 
variance spectra in Fig. Q the amplitude of the variability for the 
red wing (or the red EW) is much larger than that for the blue wing; 
hence, the basic behaviour of the total EW variability follows that 
of the red wing (the red EW). 

According to the figure, the variability of the total EW and 
that of the red EW appear to be correlated with each other. On 
other hand, the variability of the total EW and that of the blue 
wing appear to be anti-correlated for the first and the second night. 



but not on the third night. With much larger time steps Id), 
IJohns &Basri|jl995h found a similar anti-correlation behaviour of 
the Ha EWs measured using the velocity bins V ~ —150 km s^^ 
and those measured using V ~ 100 km s^^. This is consistent with 
the auto-correlation map (Fig.|6} shown earlier. 

The data points were fitted with the Fourier series (equation|4} 
for Tiniax = 2 and Jimax = 3 cases separately. As done for the 
fitting of the the subpeak positi ons in Sec tion l3.2l the period was 
constrained to be 2.7 d ( Unruh et all2004h in the fitting procedure. 
The results are also shown in Fig.0 and the corresponding Fourier 
coefficients are summarised in Table|2|along with the reduced chi- 
square (xS) values. 

For all cases (total, red and blue EWs), the lines with rimax = 

3 are better representations of the data points, as we can see from 
the figure and the values in the table. To test the validity of 
the statement above, we performed an F-test. For the total EWs 
fits, the numbers of degrees of freedom are U2 — 42 and uz = 
40 for rimax = 2 and Jimax ~ 3 (Table |2} respectively. Using 
the values of the reduced chi-square values in the table, we find 
Fx = {x^('^max = 2) - x^('^max = 3)} /x'^ = 37 and F = 
xti/xla — 8-1; hence, > F. From this test, we found that the 
inclusion of the third harmonic terms are a real improvement in the 
fitting and they should be included. Similar conclusions were found 
for the blue and the red EW variability curves. 

Interestingly, the contribution of the third harmonic terms are 
as important as the first and the second harmonic terms as one can 
see from the values of the Fourier coefficients in Table |2| This is 
consistent with the idea that variation is caused by the combination 
of the tilted-axis magnetospheric accretion model (for the second 
harmonic terms) and the presence of the cool spots of the surface at 
three different longitudes which are approximately equally spaced 
from the surface Doppler images (for the third harmonic terms), as 
briefly mentioned in Section lT2l 

Since the data sampling span of one night is very similar to 
P/3, it is very difficult to distinguish the P/3 component found in 
the Fourier analysis from a spurious detection. To overcome this 
shortcoming, a similar observation must be performed at multiple 
sites (for continuous phase coverage) for at least a few rotational 
periods. Also, authors would like to remin d reade rs that the surface 
Doppler images o f SU Aur obtained bv IPetrov et al. ( 1996) and 
iMeimessieJ il997h suffer from the same data sampling problem 
as our s. The surface Doppler images constructed bv Unruh et alj 
i2004) . on the other hand, are based on the multi-site observations 
which has a fairly continuous phase coverage and for a few peri- 
ods; however, they has significant problems with non-periodic line 
variability. 

4 MODELS 

These data, along with other spectroscopic time-series observa- 
tions, provide some strong constraints on the possible magneto- 
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Table 2. The summary of the fitting parameters in equation l^for the equivalent width (EW) data shown in Fig.Q TV and 
are the number of data points and the reduced chi-square respectively. The velocity bins used for the total, red, and blue 
EWs are -500 kms^l < V < 500 kms-^ kms"! < V < 500 kms^i, and -500 kms^^ < F < kms"! 
respectively. 
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0.25 
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Figure 7. The equivalent widths (EW) of Pa/3 using the velocity bins 
-500 kms-i <V < 500 kms-l (circles), -500 km s-^ <V < 
km (diamonds) and km < V < 500 km (squares) are 
plotted as a function time. The top panel shows the data fits using eauationl4l 
with rimax = 2, and the bottom panel shows those with rimax = 3. The 
fits for the total EW (-500 km s'^ < V < 500 km s"!), the blue EW 
(-500 km s-i < \/ < km s"!) and the red EW (0 km s'^ < V < 
500 km s~^) are shown in solid, dashed and dash-dot lines respectively. 
The rotational period, 2.7 d (Unruh et al. 2004), was kept constant in the 
fitting procedure. Tablel2lsummarises the fitting parameters. 



spheric geometry of SU Aur. Specifically, in Section l3!2l we have 
shown that the motion of the subpeak in the red absorption trough is 
related not only to the rotational period of the star, but also to half of 
the rotational period. This pheno menon is consisten t with the tilted- 
magn etic axis models (e.g. see IShuetalj[T9943 : | Johns & Basril 
Il995h as mentioned earlier. Furthermore, the analysis in Section is!?! 
suggests that the variability of the Pa/3 equivalent width is asso- 
ciated with a half and and a third of the rotational period of the 



st ar. This characteristi c is also recovered in sur face Doppler image s 
bv lPetrov et ai]^199d) ■ lMennessieJ^1997^ and|u^hetai]|2o3), 

and is indicative of the presence of cool spots on the surface at three 
different longitudes which are approximately equally separated. 

To date the observational phenomena associated with SU Aur 
have been interpreted using cartoon-like models of the circumstel- 
lar geometry and dynamics. Although useful, it is not clear that 
these necessarily simplistic descriptions provide a reasonable ex- 
planation of the changes in line-profile shape that are observed. 
However, it is also true that the near-star geometry is sufficiently 
complicated that a detailed fit to the observations is not currently 
tractable, and will probably require a combination of simultaneous 
time-series observations (spectroscopy, circular polarimetry and 
photometry) spanning a wide range of wavelength. In this sec- 
tion we adopt a modelling approach that, while falling short of 
a formal fit, provides a greater quantitative insight into the line 
profile variability than a cartoon description. Our intention is to 
develop radiative-transfer models based on the simple magneto- 
spheric geometries that have been pro posed in previous studies 
jjohns & Basrll995l:IPetrov et alll993) . We will be able to deter- 
mine whether or not these geometries are capable of reproducing 
the gross characteristics of the line profile variability, and therefore 
make a better assessment of the applicability of the models. 

4.1 The modelling code 

The three-dimension al Mo n te C a rlo radiative trans- 
fer code TORUS iHarriesI l200d iKurosawa et alj |2004 
ISvmington. Harries. & Kurosaw3 120043) is used to compute 
the Pa/3 line profiles as a function of time (rotational phase). 
First, the model computes the non-UTE populations of 14-level 
Hydrogen atoms in the magnetosphere which is funnelling the 
gas through the magnetic field lines from the inner edge of 
the accretion disc. Second, the model computes the observed 
profile of Pa/3 as a function of rotational phase. In our models, 
the Sobolev approximation (c.f. iMihalaJ Il97^ is used in the 
framework of core-plus-halo Monte Carlo radiative transfer 
method, in which the photosphere is treated separately from the 
outer atmos phere (e.g. accre t ion stre ams) of a star. Readers are 
referred to ISvmington et alj i2004ah for a detailed description 
of the spectroscopic model of hydrogen lines arising from the 
accretion streams in the magnetosphere. All of the accretion flow 
is assumed to be constant in our models i.e. there is no periodic 
pulsation of density enhancement etc. 

No rotational ve locity component is in clud ed in our models 
as it wa s neglected bv lHartmann et al] ^9^) and lSvmington et all 
<2004iJ) . The effect of the rotation on Pa/3 is expe c ted to be 
small according to the calculations of iMuzerolle et al] i200lL see 
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their fig. 8) . A slightly larger am ount of rotational broadening 



large 

than that of iMuzeroUe et al] feOOll) is expected for SU Au r since 
V sm I ~ 60 km s"^ ^Johns-Kmljllli 



Jnruh et alJl2004 while 

iMuzerolie e t al. (2001) used sin_i_« 10 km in their calcula- 
tion. Interestingly. lOliveira et^ i2000l) found the time-lagged be- 
hayiour in the variability of some optical lines, and proposed that 
this can be caused by the presence of the azimuthal component in 
the magnetosphere which might be caused by the interaction of the 
rotating magnetosphere with the circumstellar disc. Although the 
effect of the rotation on the Pa/3 line may be important for a quanti- 
tative measurement, it should not affect the conclusion drawn from 
the qualitative analysis presented in this section. The exact effect 
of the rotation on the line formation and variability will be investi- 
gated more carefully in a future paper 

4.2 Model configurations 

A star with radius i?« and mass M, is placed at the origin (O) 
of a Cartesian coordinate system {x,y,z) as shown in Fig. |8| The 
rotational axis of the star is set to be identical to the z-axis, and 
the sense of rotation is counter-clockwise when the star is viewed 
pole-on. The y-axis is perpendicular and into the page. An observer 
is placed on the z-x plane with an inclination angle i measured 
from t he z-axis. The structure of the magnetosphere is the same as 
that of lHartmmn et aljil994) i.e. the shape of the dipole (poloidal) 
magnetic field lines is described by 



r = Rn 



(8) 



where r, Rm and 6 are the radial distance from the centre of the star 
to a point (p) along a field line, the distance to the point where the 
field line intersects with the equatorial plane, and the polar angle 
measured from the magnetic axis (normally z-axis) respectively. 
The magnetic axis can be inclined by a small angle with respect 
the rotational axis. We consider the following three models: 



Model A 



Model B 



Model C 



A star is surrounded by the axi-symmetric magnetic 
field (described by equation|8} with constant accretion 
flow along the magnetic field lines, but the magnetic 
axis tilted by 10° with respect to the rotational axis. 
A star is surrounded by the axi-symmetric magnetic 
field (described by equation|8j except three thin (10°) 
gaps, where there is no magnetic field and gas, are lon- 
gitudinally placed with equal separations (i.e. 120°). 
The magnetic axis is aligned with the rotational axis. 
A combination of Model A and Model B. The magne- 
tosphere has three 10° gaps, and the magnetic axis is 
tilted by 10° with respect to the rotational axis. 



In all three models, we have adopted the following stellar pa- 
rameters: the mass M , = 2.25 M© and the radius R, = 3.6 R© 
JCohen&Kuhilll979h . The magnetospheric radius (i?m) is as- 
sumed to be comparable to the size of the co-rotation radius (Rc) of 
the star (e.g.'Pringle & Ree5^1972t lGhosh & Lambll979l : IShu et alj 



1994b : Romano va et al. 2002). Rc « 3.0 R, for i?. and M. values 
given above. The range of the magnetospheric radius is chosen to 
be Rm = 2.5 — 3.5 Rt, which is in proportion to the small/narrow 
model of lMuzerolle et all i200lh . Th e temperature structure along 
the magnetic field was adopted from lMuzerolle et alj i200il) with 
the maximum temperature (Tmax) of 8000 K (see their fig. 2). The 
stellar co ntinuum of the core star is described by a model atmo- 
sphere of lKurucJ jl979h with T^a = 5750 K and log = 4 (cgs). 

A simple geometrically thin and optically thick accretion disc 
is placed just outside of the outer edge of the magnetosphere. All 



Rotational axis 

/' Magnetic axis 



To observer 



I 



.Magnetic field line 



Figure 8. Basic model configuration. A star radius R, and mass M* is 
placed at the origin (O) of a Cai'tesian coordinate system (x,y,z). The rota- 
tional axis of the star is identical to the z-axis, and the sense of rotation is 
counter-clockwise when the star is viewed pole-on (viewed from +z direc- 
tion). The j/-axis is perpendiculai' and into the page. An observer is placed 
on the z-x plane with an inclination angle i measured from 2-axis. The 
magnetic axis is tilted by a small angle 8m with respect to the rotational 
axis. The shape of the dip ole (poloidal) magnetic field lines is described by 
r = sin-^ 6 



dip ole (poloidal) magnetic I 
e.g. lHartmann et aljfl994l) . 
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Figure 9. The Pa/3 line profiles computed at five different inclinations (i = 
90°, 85°, 70°, 60°) using the same magnetosphere as in Model A, but 
without tilting the magnetic axis i.e. Bm = 0. The subpeak in the red wing 
appeal's in emission only at the high inclination (i >~' 85°). At lower 
inclinations (i <~ 70°), the flux level of the blue side of the profile is very 
sensitive to change in inclination. On the other hand, at higher inclination 
angles (i >~ 70°), they are insensitive to the change in the inclination 
angle. 
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Table 3. Model parameters 



Model 




# of gaps 


AI 




n 




(iQ-^Moyr-i) 


A 


10 





1.15 


B 





3 


1.25 


C 


10 


3 


1.15 



the photons which encounter the disc are absorbed. With this as- 
sumption of the disc, a photon emitted from the bottom half of the 
magnetosphere can be occulted by the disc, but a photon emitted 
from the top half can not be occulted by the disc. This is reasonable 
for computing Pa/3 profiles except for cases with very high inclina- 
tion angles. If a more realistic disc (e.g. a flared disc) was used in 
the model, we would expect the reducti on of the photospheric con- 
tinuum fe.g.'C hiang & GoldreicJ 199 9*1 as well as the line emission 
due to obscuration by the outer part of the disc. Predicting the exact 
effects on the line profile shapes caused by this type of obscuration 
requires modelling of a self-consistent accretion disc, which is be- 
yond the scope of this paper. 

Since we do not have the estimate of the amount of the veiling 
around Pa/3 line, the flux contribution from the disc itself is not 
considered in our calculations. A re cent measurement of th e veiling 
at 2.2^im for SU Aur is 0.6 ± 0.3 jMuzerolle et al^2003^ ■ but the 
value should be significantly smaller at 1.3/im (Pa/3). If the veiling 
correction was taken into account, the line strength of the models 
presented in this paper would be weaker. The mass-accretion rates 
(M) for the models were chosen so that the mean observed line 
strength is approximately reproduced. 

The inclination is set to i = 80° which is higher than the 
previ ously assumed values (e.g. lMuzerolle et al .'2003':'U nruh et'al] 
120041) . Fig.|9|shows the Pa/3 line profiles computed at different in- 
clination angles using the same magnetospheric accretion structure 
as in Model A, but with no dipole offset i.e. Om = 0. At lower in- 
clinations (i <~ 70°), the flux level of the blue side of the profile 
is very sensitive to changes in inclination. This will cause too much 
variability in the blue wing if Om = 10° (as in Model A); hence it 
will be inconsistent with the observation (c.f. Figs.Qand|3}- On the 
other hand, at higher inclination angles (i >~ 70°), they are insen- 
sitive to the change in the inclination angle. Interestingly, we found 
that the higher inclination angle (i >~ 85° in Fig.|9} is needed to 
have the subpeak in the red absorption trough in emission (above 
continuum) as seen in the mean spectra of Pa/3 from the third night 
of the observation (Fig.|3)- 

In Models B and C, the locations of the gaps are at the az- 
imuth angles 105°, 225° and 345° (measured from -|-x-axis before 
tilting) at time t = or equivalently phase — 0. For the models 
with the tilted magnetic axis (Models A and C), the inclination of 
the magnetic axis with respect to the observer changes from 90° to 
70°, then to 90° as the rotational phase changes from to 0.5, then 
to 1.0. The model parameters are summarised in Table|3| 

4.3 Model spectra 

We have computed the Pa/3 spectra at 50 different rotational phases, 
and placed the results in Fig. 1101 The figure shows the mean spec- 
tra of a whole phase, the deviations from the mean spectra (quotient 
spectra) as a function of phase in the greyscale image, and the tem- 
poral deviation spectra (TVS)"^''^ from the models. 

The mean spectra of the models broadly reproduce the overall 



shape and relative intensity of the observations (Fig.0, although 
the model line profiles are rather more triangular. Little difference 
is seen in the shape of the mean spectra from three models. The 
models overestimate the levels of the blue- wing (between ~ —200 
and ~ — lOOkms^^) compared to the observation. The discrep- 
ancy could be caused by e.g. the wrong inclination angle, the wrong 
geometry of the magnetic field line used in the model, and the fact 
that the model does not include the rotational velocity component in 
the accretion streams. If a lower (e.g. 60° ) inclination angle is used, 
the profile will become less triangular. Although a lower inclination 
model will give us a better fit to the observed mean profile, the vari- 
ability in the red absorption trough will be much smaller than that 
of the observation, and the variability in the blue side will be much 
larger than that of the observation. The mean spectra from Mod- 
els A and C show the subpeak in the red absorption trough similar 
to the one in seen in the observation; however the positions of the 
subpeak in both models are 100 km s^^ smaller than that of the 
observed spectra. 

In Model A, the magnetic axis is perpendicular to the observer 
when phase = 0. Around this phase, the subpeak in the red ab- 
sorption trough become most prominent (seen as brighter pixels in 
the greyscale image). The subpeak arises from a simple geometri- 
cal effect: at high inclination the observer's line-of-sight (LOS) to- 
wards the hotspots passes through both low velocity material near 
the disc, and high velocity material near the stellar surface, leading 
to two absorption components on the red side of the line profile. 
The relative position and strength of these components changes as 
the magnetosphere becomes more face-on. The level of the subpeak 
should become lower as the phase changes from to 0.5 because 
of this geometrical/projected velocity effect. The greyscale image 
of Model A shows this effect. Overall line variability seen in the 
red absorption trough is well reproduced by Models A and C. Since 
this variability is primarily a geometrically effect, it is quite insensi- 
tive to the adopted temperature structure, and we find qualitatively 
similar behaviour for models with 6000 K < Tmax < 10000 K 
(although naturally the models differ in detail). 

The greyscale image of Model B show the periodic variation 
(with 1/3 of a phase) of the flux levels in the red absorption. The 
flux level in the red absorption increases as a longitudinal gap ap- 
proaches the line of sight of the observer since line photons suffer 
less absorption. Finally, the greyscale of Model C shows a com- 
plicated variation pattern, but it is easy to see that the pattern is 
essentially a combination of the variability patterns from Models A 
andB. 

The temporal deviation spectra of Models A and C re- 
produce approximately the same range of the line variability 
(-400kms"\ 500 km s"^) as seen in the observations (Fig.0. 
Model B does not fit the (TVS)"^/^ of the absorption at all. In 
the (TVS)^^'^^ of Models A and C, rather large discrepancies are 
seen in the amount of the variability around the line centre (slightly 
on the blue side). While the peak of the variability in the red ab- 
sorption occurs at ~ 150 km s^^ for Model A and C, it occurs at 
~ 275 km s~^ in the observation. 

A possible way to shift the peak of the (TVS)^^''^ in the red 
wing (equivalently the position of the subpeak in the mean spectra) 
to a higher velocity bin is to have a larger magnetospheric radius 
(i?m), in which the accretion gas has higher infalling velocity. Al- 
though this means the magnetic accretion streams have to extend 
beyond the co-rotation radius of the star, this w ould be still con- 
sistent with the magneto-hydrodynamic model of iRomanova et all 
i2002h who found a torque-less accretion is possible when i?m ~ 
1.5i?c where i?c is the co-rotation radius of a star. 
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Figure 10. The summary of the Pa/9 spectra computed for Model A (left). Model B (centre) and Model C (right). For each model, spectra were computed at 
50 different rotational phases. In the bottom panels, the mean spectra of all rotational phases are shown. In the middle panels, the quotient spectra (divided 
by the mean spectrum) are shown as greyscale images with increasing rotational phases in vertical direction. The greyscale image is scaled from 1.1 (white) 
to 0.9 (black). The temporal deviation spectra (TVS)~^/^ are shown in the top panels. When computing the temporal deviation spectra, the signal-to-noise 
ratio of 90 (in continuum) was used to match that of the observation (Section|2|. See Tablel3lfor the model parameters adopted. 



4.4 Model auto-correlation maps 

The auto-correlation maps of Pa/3 line (with no time lag) were com- 
puted using the spectra of Models A, B and C, and the results were 
placed in Fig. II II along with the auto-correlation map shown in 
Fig.|6|for a comparison. 

In the map for Model A, the changes in the flux 
levels in the velocity range (—200 km s^^, km s^'^) and 
(0 km s~^, 400kms~^) are strongly correlated with themselves 
while they are anti-correlated with each other. In other words, as 
the flux level in (—200 km s~^, km s~^) increases the flux level 
in (0 km s~^, 400 km s"'^) decreases, and vice versa. 

The blocks of correlation and anti-correlation (seen 
as black and white squares in the map of Model A) are 
less pronounced for Model B. The third quadrant block 



(-300 kms" 



0kms"^)x(-300 kms" 



km s"^) of 



Model B is more complicated than that of Model A because the 
accretion streams located on the far (back) side of the star (as seen 
by the observer) will be seen through the gap as it approaches the 
near (front) side of the star. The map for Model C is essentially 
the same as that of Model A since the overall variability caused by 
the presence of the gaps is relatively small compared to that by the 
precession of the magnetic axis. 

The main differences between the auto-correlation map of 
the observation and the models are: 1. the absence of the anti- 
correlation between the velocity ranges (Okms"\ 100 kms"^) 
and (200 km s"^, 400 km s~^) in the models and 2. For Models A 
and C, the two dark stripes (anti-correlation bands) crossing per- 
pendicular to each other are about two times wider than that of 
the observation, and they appear as two large dark blocks and two 
smaller blocks. The total phase coverage of the observation is ap- 
proximately 1 /3 of the whole phase; therefore, this may also con- 
tribute to the difference seen in the maps from the observation and 
the models. 




Figure 11. Auto-correlation maps of Pa/3 from the observation (upper- 
left) and from the radiative transfer models: Models A (upper-right), B 
(lower-left) and C (lower-right). The velocity bins of the continuum are 
excluded from the model auto-correlation maps, and their values are set 
to (seen as gray outer frames). No anti-correlation between the velocity 
ranges (0 km s^^, 100 km s~^) and (200 km s~^, 400 km s~^) is seen 
in the models unlike the observation. The position of the two dark stripes 
(anti-correlation bands) crossing perpendicular to each other is on the blue 
side (—125 km s^^ and —125 km s^^) for Models A and C, but it is on 
the red side (100 km and 100 km s~^), in the observation. 
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4.5 Motion of the subpeak in tlie red absorption 

In Section lT^ we found that the motion of the subpeak in the red 
absorption trough was possibly related with the precession of mag- 
netic axis around the rotational axis. In Fig. 1121 a subset (25 out 
of 50) of Pa/3 line spectra from Model A is shown. Only the red 
side of the line profiles are shown to emphasise the motion of the 
subpeak. At phase = 0, the inclination angle of the magnetic axis 
with respect to an observer is 90°, then it gradually decreases to 
70° as the phase reaches to 0.5. The inclination of the magnetic 
axis increases again to 90° as the phase reaches to 1. 

The velocity position of the subpeak moves from V ~ 
120 km s"^ (at phase = 0) to 1/ f« 200 km s"^ (at phase = 
~ 0.25), and then to V ^ 130 km s"^ (at phase = 0.5). The 
motion is symmetric about phase = 0.5. As the phase approaches 
0.25, the main subpeak feature becomes less pronounced, and it 
smoothly shifts to a smaller subpeak seen at V ~ 130 km s~^ (in- 
dicated by an arrow in the figure) as the phase approaches 0.5. 

Fig. 1131 shows the positions of the subpeak, presented in 
FigEI 3S 3 function of rotational phase. It clearly shows the 
presence of the second harmonic components as seen in the ob- 
servation (Fig. |4}. At phase around 0.25 and 0.75, the velocity 
reaches the maximum value {V ~ 200 km s^^) which is about 1.4 
times smaller than the maximum velocity seen in the observation 
{V ~ 280 km s^'^). At phase = and 1, the velocity reaches the 
minimum value (V ~ 120 km s~^) which is comparable to that of 
the observation (V « 140 km s^^). 

In the time-series spectra of Model B, no clear subpeak is seen 
in the red absorption trough therefore it is not shown here. Model C 
shows a similar variability pattern of the subpeak as in Model A. 
The gaps in the magnetosphere in Model B and C, are unlikely the 
causes of the subpeak motion. From this analysis, we found that 
the Models A and C reproduce qualitatively similar variability in 
the red absorption as seen in the observation. 



4.6 Line equivalent width variability 

Based on the model spectra computed in section l43l the EW of 
Pa/? were computed using the velocity bins from —500 km s^^ to 
500 km s^^, and plotted as a function of time (Fig. ll4>. The model 
phase was multiplied by the period of SU Aur (2.7 d JUnruh et all 
l2004h . and translated by -1-9.1 d to match approximately the posi- 
tions of the peaks in the observed EW variability curve. 

According to Fig. ll4l the ranges of the EW variation predicted 
by Models A, B and C are about 0.7, 0.5 and 0.9 A respectively. 
On the other hand, the observation (Fig.|7} shows the ranges of 0.7, 
0.5 and 1.2 A for the first, second and third nights. The amplitudes 
of EW changes predicated by the models are within the range of the 
observational uncertainty (~ 0.3 A). A direct comparison between 
the absolute values of EW would be misleading since our model is 
not a formal fit to the data. Nonetheless, we note that our model 
mean EW (—1.0 A) is comparable to the mean EW of the obser- 
vation (—0.7 A), although of course the observations do not fully 
sample a complete rotational period. 

The overall characteristics of the EW variability seen in the 
observation is well reproduced by Models A and C, but Model C 
with three gaps is preferred because it shows both P/2 and P/3 
variation found in the observation (section|T4). Model B is clearly 
inconsistent with the data. A better match of the observation with 
Model C can be achieved by carefully choosing the longitudinal 
position of the gaps with respect to the tilt of the magnetic axis. The 
amplitude of the EW variability can be increased by using wider 
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Figure 12. The red wings (V > km s~^) of tlie Pa/3 spectra from 
Model A at 25 different rotational phases (0-1). The velocity position of 
the subpeak moves from V ~ 120kms~^ (at phase = 0) to V Ri 
200kms"i (at phase = 0.25), and then to V ^ 130 kmg-^ (at 
phase = 0.5). The motion is symmetric about phase = 0.5. As the 
phase approaches 0.25, the main subpeak feature becomes less pronounced, 
and it smoothly shifts to a smaller subpeak (indicated by an arrow) seen at 
V ~ 130 km s~^ when the phase is around 0.5. Also see Fig. 1131 




Rotational Phase 



Figure 13. The positions of the subpeak in the red wing of the Pa/3 time- 
series spectra (Model A) shown in Fig ll2l are plotted as a function of rota- 
tional phase (0-1). The figure clearly shows the presence of the second har- 
monic components as seen in the observation (Fig.0. At phase around 0.25 
and 0.75, the velocity reaches the maximum value {V ^ 200kms-l) 
which is about 1.4 times smaller than the maximum velocity seen in the 
observation (V ~ 280 km s~^). At phase = and 1, the velocity reaches 
the minimum value {V ~ 120 km s~^) which is comparable to that of the 
observation (V ~ 140 km s~^). 
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Figure 14. The equivalent widths (EWs) of Pa/3 as a function of time from 
Model A (top), Model B (middle), and Model C (bottom). The mo del phase 
was multiplied by the period of SU Aur (2.7 d JUmuh et all2004 . and then 
translated by +9.1 d to match the peaks of the observed EW. Models A and 
C reproduced the overall characteristics of the EW variability seen in the 
observation (see Fig.0, but Model B does not. 

gaps in the magnetosphere, and by increasing the tilt angle of the 
magnetic axis. 

Interestingly , the recent light curve calculation of 
iRomanova et alj i2004l) . who considered the flux variation 
caused by the hot spots on the surface of the rotating stars in 
their magneto-hydrodynamic simulations, showed the presence of 
second harmonic (two peaks in one rotational period) when the 
inclination angle (i > 75°). This is consistent with the EW curve 
of our Model A. 



5 CONCLUSIONS 

We have presented the 503 time-series echelle spectra of Pa/? line 
of SU Aur (Fig. 0. From the temporal variance spectra analysis 
(Section |3.1> , we found that the line displayed relatively strong 
variability within the velocity range (100kms~^, 420kms~^) in 
the red absorption trough, and weak variability within the veloc- 
ity range ( — 200 km s ~ ^ , km s ~ ^ ) in the blue wing while the flux 
level at the line centre remained almost constant. A hint of blue- 
shifted absorption at about — 150kms^^ was found on the sec- 
ond night. Interestingly, a similar blue-shifted absorption feature 
at — 150kms^^ are seen in Ha and H/3 profiles, which modu- 
lates with the rotati onal period of the star iGiampapa et alj[l993l : 
I.Tohns&Basrll99.5ft . 

We found the presence of a subpeak in the red absorption 
trough in each spectra (Fig. |3}. and its position changes in time 
(Fig. |4}. The variability seen in the red-wing seems to be caused 
mainly by the motion of this subpeak. The position of the peak was 
fitted well with the Fourier series summed up to the second har- 
monic term (rimax = 2 in equation|4}, indicating that the variabil- 
ity is associated not only with the rotational period but also a half 



of the rotational period. This is co nsistent with the tilted-axis mag- 
netos pheric accretion models (e.g. lShu et all 1994^ : 1 Johns & Basril 
Il995h . To exclude the possibility that the large amplitude of the 
change in the subpeak position seen on the second night was caused 
by a single episodic event, the object must be observed for a few 
rotational periods. 

The auto-correlation map (Fig.|5} showed the profile variabil- 
ity in the velocity range (— 200 kms^^, km s^^) weakly corre- 
lates with that for (200 kms~^ , 400 kms~^). The variability in 
those two velocity ranges weakly anti-correlates with the variabil- 
ity in (0 kms"'^, 100 kms~^). The map is v ery similar to that see n 
in the H/3 auto-correlation function shown bv lOliveira et aljfeood) . 

The modulations of the line equivalent widths (EW) (Fig.0 
are fitted well with the Fourier series (equation|4} summed up to the 
third harmonic terms (rimax = 3 in equation|4}, indicating that the 
modulations might be associated not only with the rotational period 
of SU Aur and a half of the period, but also a third of the period. 
Because our data sampling period and the P/3 component of the 
Fourier analysis are very similar to each other, it is possible that 
the detection of the P/3 contribution may be spurious. A definitive 
observation requires a multi-site campaign to continuously monitor 
the object for at east two rotational periods. 

From the comparison of the observation with the radiative 
transfer models (Models A, B and C in Section|4j, we have learnt 
the foUowings: 

1. The inclination angle (with respect to the magnetic axis) 
must be rather high (~ 85° ) to have the subpeak (in the red wing) in 
emission as seen in the mean observed profile from the third night 
(see Figs.|3|and|9j. If a more realistic accretion disc was used in our 
model, instead of the geometrically thin and opaque disc, this angle 
would most likely change. At the lower inclination angels (i <~ 
70°), the models will overestimate the variability in the blue wing 
when compared to the observation. Since the profile shape highly 
depends on the geometry of the magnetosphere, these conclusions 
would change if the geometry other than the one in Hartmann et al] 
il994 is used. 

2. The tilted magnetic axis models (Models A and C) repro- 
duce the temporal deviation spectra (TVS) ^ ^''^ which have similar 
ranges and shapes as seen in the observation (Fig. llOt . On the other 
hand, the mean spectra of the models (A, B and C) are more asym- 
metric and triangular than the observed profile. The discrepancies 
may arise from several simplifying assumptions in the model, in 
particular the adoption of Sobolev line transfer (which is a poor 
approximation at the line centre) and the neglection of turbulent 
broadening. 

3. The motion of the subpeak in the red absorption trough seen 
in the observation (Figs.|3|and|4} was qualitatively reproduced by 
the tilted magnetic axis models (Models A and C). See Fig. 1121 The 
models reproduced the motion of the subpeak with the second har- 
monic components (associated with P/2 where P is the rotational 
period), as seen in the observation. 

4. The tilted magnetic axis models with the gaps in the magne- 
tosphere (Models C) predicts the line equivalent width (EW) vari- 
ability patterns very similar to that seen in the observation (Figs.Q 
and ll4t . The model shows the P/2 and P/3 variability of the EWs 
as seen in the observation. 

The models presented here represent the next step in quanti- 
fying the magnetospheric geometry of T Tauri stars, moving from 
purely qualitative descriptions of the circumstellar environment to 
a more quantitative approach. The complex nature of the problem, 
with inflowing and outflowing material, non-axisymmetric (and 
possibly non-stationary) accretion, render a formal fit to the data 
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extremely difficult. Further progress will only be possible with an 
extensive, coordinated observational programme, spanning a vari- 
ety of techniques and wavebands. Zeeman Doppler Imaging (ZDI; 
[Semeletal. 1 991 of the surface magnetic field topology would 
provide strong constraints on the magnetospheric geometry, par- 
ticularly when combined with photometric time-series. Simulta- 
neous Hq and Pa/3 spectroscopy would also be extremely useful, 
since Ha is appears to be a good tracer of outflow emission, while 
we have shown that Pa/3 probes the accretion streams. Time-series 
Hq linear spectropolarimetry may also provide useful new insights 
jVink et alliooilvink et all2004 . 
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